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A model is developed to explain the temperature dependence of the group velocity as observed 
in the experiments of Hau et al (Nature 397, 594 (1999)). The group velocity is quite sensitive to 
the change in the spatial density. The inhomogeneity in the density and its temperature dependence 
are primarily responsible for the observed behavior. 
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I. INTRODUCTION 



The phenomenon of Bose-Einstein condensation in atomic gases |l|J^ lends itself to the study of many fundamental 
effects. Among them, one aspect presently being investigated both theoretically and experimentally is the interaction 
of light with atoms in the quantum degeneracy regime 1^. In this context, the propagation of light inside a cold 
gas is still an open problem. Because of the optical density, it is well known that the transmission of resonant light 
through a condensate is almost zero However, electromagnetically induced transparency (EIT) Q was found 
to allow the propagation of light by means of quantum coherence between different internal atomic levels |^,^. In 
this context, Hau et al discovered a remarkable property of pulse propagation in a Bose condensate. These authors 
demonstrated the slowing down of the group velocity of the pulse to 17 meters/sec Furthermore, they have 
shown a definite dependence of the group velocity on the temperature of the ultra-cold sample. One would like to 
understand the observed temperature dependence from first principles. For this purpose, it is necessary to extend 
the standard theory of EIT to a cold gas at finite temperature. However, a theoretical description of this problem 
is rather complex. Complexities arise when one attempts a systematic treatment of interactions, finite temperature 
effects and dynamics. Most studies treat these aspects as disjoint: interactions are included in the zero temperature 
case to study the kinematical aspects whereas some dynamical aspects are studied using only the excitations 

within the electronic ground state and finite temperature effects are usually studied for noninteracting bosons 
|p|, P^|JTl]] . A complete theory should study all these aspects together. However, a complete theory of the interaction 
of light and interacting particles is still unavailable, and a full numerical treatment is a rather hard task. Here, we 
present approximate but plausible arguments to explain the experimental observations in |^ . The simplicity of our 
model allows for an analytical expression for the group velocity in the following cases: atoms confined in a box and 
by a harmonic potential. We obtain results which reproduce the ones in for T > T^. In particular, the treatment 
brings out the factors playing key roles in the phenomenon. Here, we show that the variation of spatial density of 
atoms with temperature is the major factor responsible for the temperature dependence of the group velocity. 
The paper is organized as follows: In Section II the model is introduced. In Section HI we derive the group velocity 
of a pulse propagating in an ideal gas confined inside a box, extend the calculation to the case of an ideal gas in 
a harmonic oscillator potential, and present and discuss the results in relation to the experiment of Hau et al. In 
Section IV we present estimates for the group velocity in the interacting case and in the limit of zero temperature. 



II. THE MODEL 



In this Section we introduce the model used throughout this article. Here, we write the Maxwell-Bloch equations 
which describe the dynamics of the system consisting of light field and atoms. We derive the linear response of the 
medium to a weak probe field, taking into account the quantum statistics of the atoms. The group velocity is then 
defined in the standard manner |12(|. 



1. Maxwell-Bloch Equations 



We consider a gas of N non-interacting bosons. The relevant internal structure corresponds to a three-level atom, 
with internal levels \g) (stable state), |r), (metastable state) and |e) excited state, whose energies are Wg, ujr and We, 
respectively (see Fig. 1). The radiative decay rate of the excited state is 7 = 7g -1-7^, with 7g(7r) the rate of decay on 
the transition |e) |ff)(|e) \r)). Laser light with frequency wig and wave vector kg drives the transition \g) |e), 
whereas the transition \r) — > |e) is driven by a field of frequency uJir and wave vector kr- The dynamics of the whole 
system is given by the Maxwell equation for the electric field vector E 
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and by the optical Bloch equations for the density matrix equations of the iV-atom gas. For non-interacting atoms 
it suffices to consider the equations for the one-atom density matrix p, projected on the basis e)} with j = r, g,e 
and |e) the eigenvector of the mechanical motion of one atom at the energy e. They have the form: 
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^Prr{e, e) = ~i{e - e')Prr(e, e') (3) 
^Pee(e, e') = - e' - i7)Pee(e, e') (4) 
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where p.y(e,e') = (i,e|p|j,e') = r,e,5), C^^^^, = (e|exp(ikj • r)|e'), = peje-'^'^*, p^g = p^^e-'^'^'s-'^"-)* and 
Pij — (Pji)* for i 7^ j'. Here, c is the speed of hght, P is the polarization of the medium, A^ ~ lu^ — Wj — ujij {j ~ g, r) 
is the detuning. Rabi coupUngs are given hy g — |deg • E|/?i and Q = |der • E|/?i, where dgj is the dipole moment 
of the transition |e) — > Finally, p^g(e,e') describes the density matrix after a spontaneous emission event on the 
transition |e) 
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= E E / d%|dje-a'(k)p(e|e*J-|ei)pee(ei,e2)(e2|e-'^^-n6') , (8) 
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where a^(k) and a^(k) form a set of polarization vectors orthogonal to k, and dje — dje/|dje|. 

In deriving the above equations, we made the rotating wave approximation and transformed to a reference frame 
rotating at the optical frequency of the laser. Furthermore, in Eqs. (||)-(|7|) we have introduced the loss-rates F^j, 
which take into account the effects of other mechanisms of decoherence. In the ideal case Fge = F^e = 7/2, whereas 

Tgr = 0. 

2. Susceptibility and Group velocity 
In an isotropic medium, the linear susceptibility x is defined by the expression 

p{t)= f dt'x{t-m{t') . (9) 

Assuming that two light fields are propagating through the medium along the z-direction, we can write the electric 
field E{y, t) and the atomic polarization P(r, t) as E{y, t) — J2j=g A^oi^^ exp{+ikjZ — iujijt) +c.c.]/2 and F(r, t) — 
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J2j=g rl-^o exp(-|-ifcjZ — iuJijt) + C.C.J/2 where Eq, Pq are respectively the slowly-varying envelopes of the electric 
field and atomic polarization at frequency coij. The Fourier transform of Eq. (^) gives Poij-o) = EQ{Lu)x{io). From 
the relation P{t) = Tr{p^(<)d} for the polarization with d the atomic dipole moment operator, we find that the 
macroscopic polarization of the medium at the position r is 

P^ir,u;i,)/2 = {r\p%dg,e-"^^^\r) , (10) 
where is the A^-atom density matrix, that has the following form in the energy representation. 

4 = E4(^'0k)(e'l . (11) 

The N-atom optical coherence density matrix has to be obtained from the solution of Eqs. (||)-(@) subject to the 
initial condition: 

p^(0) = ^7V(e)|ff,e)(ff,c|, (12) 

e 

where iV(e) is the number of atoms in the ground state with energy e. The Eqs. (|l)-(0) are to be solved to first order 
in the field Eq and to all orders in the field Eq. In this work we are interested in the steady state of the atoms with 
the field, which is assumed to be reached on a time-scale much shorter than the thermalization time-scale of the gas. 
On the basis of this hypothesis, we assume that the initial condition ( p^ ) and Eqs. (H)-(0) determine the steady state 
solution. Note that the coefficients Cg_e' determine the one-atom energy states involved in the transition induced by 
the laser field. For free bosons these coefficients have a simple form, as it can be seen in section III. A. For harmonic 
oscillator potentials the coefficients C^.e' are given in terms of Laguerre polynomials, where the number of vibrational 
states which are coupled depends on the ratio between the recoil frequency over the trap frequency. In section III.B 
we use an approximate treatment for this case. 

Once the susceptibility is known, the dispersion relation of light in the medium is given |l^ and we can evaluate the 
group velocity, defined as Vg = ■§^\(^=ujcg- In the limit Nx <C 1 the group velocity has the form: 

- aZ' ' (13) 
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where x' = R'C(x). 



III. EVALUATION OF THE GROUP VELOCITY 



In this section we derive an analytical expression for the group velocity of a laser pulse propagating through an 
ideal gas of ultracold atoms. We investigate two cases: atoms in a box and atoms confined by a harmonic oscillator 
potential. Finally, we do the numerical calculations for the case of a gas of sodium atoms and discuss the results in 
relation with the experimental data of M . 



A. Group velocity in a gas of free non-interacting bosons 

We consider a gas of N bosons in a box of volume V. In this case the atomic wave vector eigenstates |k) are also 
energy eigenstates with eigenvalues e — where m is the atomic mass. So, we project Eqs. (|])-(0) on the motional 
basis {|k)}. Then the coefficients C^^, appearing in the density matrix equations have the form C;^^, = Cj^ ~ (^k.k+kj ■ 
We substitute these values into Eqs. (H)-(0) and solve the equations in the steady state limit. To first order in g/il 
and g/F, and assuming that at t = the gas is in thermal equilibrium, the steady-state optical coherence peg is found 
to be 

Pe.(k-k„k)--^ 2.(r,.-fz(A,-A.)) 



{Tge + iAg) {Tgr + l{Ag - A,)) 

(14) 
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where Ao is the detuning defined as 



A,=A," + ^^ , (15) 



with A° = A^+Ufi, and ujr is the recoil frequency defined as luh — hk"^ /2m. Using ( JI^ in (p^), we find the expression 
for the susceptibihty 

X,eK)-X°E^W-^ A„ . . ^ OV4 ■ (16) 

The sum in Eq. is over all the motional states weighted by their statistical occupation A^(k), 

/ ^ exp ypfi k^/2m) — 1 

where / is the fugacity, (3 = 1/KbT and T is the temperature, and k — |k|. Here, x° is the one-atom susceptibility, 
defined as 

_ |dgeP _ , . 

^ - 1^ = 32^ ' 

where A is the optical wavelength of the transition g ^ e, X — 27rc/a;ge, and Tge — 7/2. In the following, we assume 
that |kg — kj.| <C kg, kr, Tgr- Therefore, the dependence on k in the denominator of Eq. ( p^ ) is mainly due to the first 
term A^, and we may rewrite Eq. jl^ ) as 

where A;^ = k • i and ^ is a complex number independent of fc; 

A° y4 

^^"^i^'^V,4^,. + .(A0-A0)) • (2°) 



After evaluating the susceptibility as given by Eq. (g^, the group velocity can be found using Eq. (|13D. 

In the following, we investigate the behaviour of the cloud close to the critical point, dividing our investigation into 

two regimes: above and below the critical temperature Tc- 



1. Above the critical temperature 



Above the critical temperature and in the limit of large volumes, one can replace the sum in Eq. ( jl^ ) with an 
integral in three-dimensions jl^ . The expression to evaluate is now 



(2^)3 
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where for convenience we have chosen to integrate in the cartesian coordinates. We write iV(k) as 

oo 

N{k) = XI /' {"Wfi^k^/2m) , 

and use it in Eq. (pl^) to obtain 
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This can be written in terms of the standard functions 

3/2 OO ,1 



where 



2KbT kg 



and where the function w is defined as 

w{x) = exp(— x^) (erf(ia;) + 1) 
Given the critical temperature for an ideal Bose gas 
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where n(= N/V) is the density of atoms, we can rewrite Eq. ( p4[ ) as 

f ^ ■ ^ 1 V- /' ( 
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2. Below the critical temperature 
For T < Tc, we use the expression for the ground state population in the thermodynamic limit ||14|, to obtain 
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Xge{^^lg)+\f=l - —n 
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where the second term on the RHS describes the contribution of the condensed phase. 
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3. Regime of parameters and approximations 

For currently studied optical transitions the argument of the w function in Eq. (Eq) is y — \/lC,l A S> 1 for any 
value of I > 1. Therefore, the asymptotic expansion of the w function can be applied |16[|: 
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We substitute this expansion into Eq. (^8|), and obtain for T > 

.T A 



Tc y/^Ac g3/2(l)C ^ 
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where we have used the relation 53/2(/)/33/2(l) = Tc/T for T > Tc- For T < Tc the susceptibihty will now have the 
form 



Xgeii^lg)- = 



l+( n'^'. 95/2(1) 



53/2(1) 



(34) 



Using Eqs. (|33| ) and (|34[)in the formula (13) we find the group velocity. Note that the dependence on the temperature 
comes in at higher order in the expansion A/(^. Clearly a significant temperature dependence for a free gas can come 
only for narrow optical transitions. 



B. Group velocity in a gas of trapped non-interacting bosons 

Let us now consider a cloud of atoms trapped in a three-dimensional harmonic potential with cylindrical symmetry, 
so that the one-atom Hamiltonian describing the mechanical motion has the form 

H^^ + Vir), (35) 

where V(r) is the harmonic oscillator potential in cylindrical coordinates 

Vir)^^m{iyy + iy^,z^) , (36) 

with Vr, trap frequencies in the radial and axial directions, respectively. In order to evaluate the susceptibility in 
the steady state, we solve Eqs. (|2])-(^ in the semiclassical limit for the atomic motion, and we sum over the states 
using the semiclassical statistical distribution [p|JT^]. This limit is valid when treating the non-condensed fraction of 
atoms for temperatures T fulfilling the condition KbT ^ fiv^ and under the condition F, A 3> i'. The hypothesis is 
justified in the range of parameters of and simplifies considerably the treatment, allowing for an analytical solution 
of the group velocity. Then the coeffients C^^^' simplify to their semiclassical values {Cc,e') ~ ^p,p'-(-fik^r,r', where p, 
r are now the classical canonical coordinates of a harmonic oscillator with energy E — p^/2m -I- V{y). In this limit, 
the optical coherence Peg appearing in Eq. ( p^ ) has the form: 

Peg(r,p) = -^^^ . . ^ ' (37) 

r,e 'T,,(V,r+t(A,-Ar)) 

with Aj defined in Eq. (Iq), and the susceptibility is given by the expression: 



d^p 
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where the semiclassical statistical distribution is 



A^(r,p) = 1 . (39) 

/-iexp(;3(i^ + y(r))-l 
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When considering the condensate contribution to the optical susceptibihty, one should evaluate Peg{€,e') and sum 
over the final states with energy e. However, in the regime lvr/v 3> 1 we may apply the semiclassical approximation 
to the final states. The final semiclassical energy is the recoil energy, and we can write the optical coherence for the 
condensate contribution as 

r^e *r,e(r,.+»(AO-A?)) 
The ground state occupation is given in the thermodynamic limit by 

(41) 

where Tc is the critical temperature of the trapped gas, KbTc = h{i'zi^fy^^{N/g3{l)y^^. 

Contrary to the case of free bosons, the group velocity is not directly given by the formula (|3|), because of the spatial 
variation of the atomic density and therefore of the susceptibility. Here, we evaluate the group velocity using a method 
equivalent to the experimental one of i.e. we estimate the size D of the cloud and calculate the delay At of a 
pulse propagating across a selected region of a cold gas with respect to a pulse propagating in the vacuum. The group 
velocity is then given by the ratio of the size over the delay Vg^^ = D/At. Assuming that the light is propagating 
along the z-axis and cuts a cylinder inside of the volume with section S and centered on the z-axis of the cloud, we 
write the delay {At) as the average over the section S of all the delays At{r) of pulses propagating at distance r from 
the axis of the cloud 

(At) - ^ / dr2nrAt{r) , (42) 
where R is the radius of the illuminated circular section S of the cloud and At{r) is defined as 

rL{r) 

At{r) = / dz[vg{r,z)]-^ , (43) 

where Vg is defined in Eq. (p^), and L{r) is half the length of the path along the cloud. Note that R will, in 
principle, depend on the size of the incoming Gaussian beam. However, in the experiment of Q, i? is the radius of a 
pinhole set before the measuring apparatus. The delay time is experimentally obtained by measuring the difference 
between the delay time of the pulse propagating across the cloud and the one of a pulse propagating in the vacuum. 
Assuming that L(r) = L is the distance between a slit before the cloud and the photomultiplier, the final delay will 
be At = {At) — L/c. In the following, we evaluate the group velocity as a function of the temperature above and 
below the critical temperature. 



1. Above the critical temperature 



The integral over the momenta in Eq. (|38| ) can be evaluated along the lines of the procedure outlined in Eqs. 
(|l|)-(||). One finds 

Xge{r,UJig)+ (44) 

Again, the considerations on the w function made in the free case are applicable, and using its asymptotic expansion 
[Eq.(|32l)], one gets 

Xge(r,^ig)+ = -— ( ^_^2 ) (45) 
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For L ^ Dz{T), where Dz{T) is the axial thermal size of the cloud, we can replace L by cxd in the integral (|4^). 
Therefore, the delay of a beam propagating along the z-axis, and entering the cloud at a distance r from the cloud 
axis is 



At{r)+ = -2n 



1 



cdAe 2^n\, 
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1=1 



(46) 



Above the critical temperature the total delay (At)+ is thus 

LodCx'^iKBTf 2 



(At). 



-2n 



c dA e n^v.v^r 



1 



1 - e 



(47) 



Here, we take the size of the cloud to be the variance of the thermal distribution along the z-axis, and thus D^, ~ 
yj2KBT jmv^. The group velocity above the critical temperature is then Dz/{At)j^. Note that, in the limit R <C 

Dr{T), where Dr[T) = ^2KBT/mv^ is the radial thermal size of the cloud, the exponential appearing in Eq. ( [47D can 
be expanded to yield (A<)+ w l/T. Since D^iT) cx -v/T, the group velocity depends on the temperature as Vg (x T'^/^. 
The same behaviour can be found when considering the other limiting case, i.e. R ~ Dr{T) / ^/2, corresponding to 
averaging over the whole cloud. Then, the delay time has the form: 



(At>- 



-ZTT 



^4 ( 1 



c dA C2 h^iy,„^ nR^ ^ P 
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= -27r -—N 
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e 53(1) 

e 53(1) 



(48) 



where we have used the definition of critical temperature and the relation 53(l)/(?3(/) = (T/Tc)^ for T > Tc- From 
Eq. (|4^) one sees that the dependence of the delay time on the temperature of the sample appears principally in 
the spatial-average term 1/ R?, which is proportional to 1/T. Thus, the main dependence of the group velocity on 
the temperature comes in through the volume of the cloud, since Vg = D,{T)l{At)+ cx D,{T)R^ cx T^/^, and the 
variation of the group velocity with temperature is mainly due to the change of volume of the cloud. 



2. Below the critical temperature 



For the ground state of the harmonic potential, we use the expression for the ground state population in the 
thermodynamic limit ||^ and the optical coherence as given in Eq. ( pO[ ) to obtain 



Xse(r,W/g)_ =Xge{r,^^lg)+\f= 



(49) 



where N is the total number of particles, and (j){v) is the harmonic oscillator ground state wavefunction. The delay 
(At)_ is given by 

(At)_ = (At)+|y=i + (Ai)c, (50) 
where {At)c is the contribution to the total delay given by the pulses which cross the condensate 



(51) 



Here Fc is the average of the ground state wavefunction, which, according to Eqs. (Bq), (tta), (W3) is 
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where aoj = ^Jfi/mvj is the size of the ground state of the harmonic oscillator in the j direction (j = r,z). In order 
to evaluate the size of the cloud, we consider that for T < Tc, a fraction {T/TcY of the atoms is outside of the 
condensate, whereas a fraction 1 — (T/Tc)^ is in the condensate. Applying the semiclassical approximation to the 
non-condensate part, we have that 



(z2) = {T/T,f{z^)MC + [1 - {T/T,f]{z^)c, 
which leads us to defining the size of the cloud to be: 



1 



-'Oz 



Dividing by (At)_ we find the group velocity below the critical temperature. 



(53) 



(54) 



C. Numerical Results 



In Fig. 2 we plot the group velocity of a gas of sodium atoms in a box (dashed line) and in a harmonic oscillator 
(dotted line) as a function of temperature T, scaled according to the critical temperature of each case. Density of 
atoms, number of atoms and trap frequencies have been taken from the data of [Q. The experimental results of Q are 
seen to be broadly in agreement with the harmonic oscillator case. The inhomogeneous spatial density of the atoms 
and its variation with temperature is the key to the understanding of the experimental data. The curve representing 
the case of free atoms shows that the temperature dependence entering into Eq. (|2^) as a higher order correction has 
a negligible effect on the considered scale, and cannot be interpreted as the cause of the behaviour observed in . 
In Fig. 3, we compare the group velocity for two different values of the Rabi frequency coupling |r) to |e). The 
behaviour for T > Tc is similar to the corresponding one measured in |^] . The curves we obtain are however steeper, and 
this can be explained by considering the approximations made in our treatment. In our calculations we have assumed 
the same number of atoms at every temperature. However, in the experiments lower temperatures are achieved 
by means of evaporative cooling. This implies that the points of the experimental curve at higher temperatures 
correspond to larger numbers of atoms, and correspondingly to larger spatial density (for ideal gases). This leads 
to a smoother gradient of the group velocity versus the temperature than in our case. On the other hand, as the 
temperature decreases, the effect of the interactions gets stronger causing, among other effects, a lower density of the 
atoms than in the non-interacting case. Hence, one would expect a group velocity value larger than the evaluated 
one. Albeit these considerations, the evaluated curve reproduces the experimental one above the critical temperature 
with some agreement, showing that the ideal gas model provides a qualitative description of the phenomenon. 
Below the critical temperature the discrepancy between the experimental data and our theoretical predictions is rather 
dramatic. This is not surprising since the size of the condensate is strongly affected by the effect of the interactions. 
Already Ketterle and co-workers have reported that the cloud size is much larger in the interacting system compared 
to the size of the harmonic oscillator ground-state wave function p8[ |. Therefore, our evaluation can be expected to 
lead to smaller values of the group velocity than the experimental records. In order to illustrate this point, in the 
following section, we estimate the group velocity at T = by comparing the ideal case with the Thomas-Fcrmi case. 
Finally, we discuss the measure of the group velocity in the two limiting cases for a section S with radius R <^ (T) 
and with radius R « Dr{T). This is illustrated in Fig. 4, where the same dependence of the group velocity on the 
temperature is evident. The orders of magnitude of the pairs of curves corresponding to the same set of parameters 
are comparable, showing that the behaviour observed in [Q originates mainly from a change in the "average" spatial 
density of the gas with temperature. 



IV. GROUP VELOCITY FOR AN INTERACTING BOSE GAS 



In this section, we compare the group velocity value at 
the particles as non-interacting, and the interacting case, 
estimate the group velocity using the set of parameters of 



r = in the two limits: the ideal one, where we consider 
which we treat in the Thomas-Fermi approximation. We 
the experiment and the formula (1) of 0: 
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(55) 



where n is the density. Therefore, we need to evaluate the group velocity at T = by substituting into Eq. ( |55| ) an 
estimate of the spatial density, which we calculate here as the ratio of the total number of atoms over the volume of 
the cloud. This evaluation, which corresponds to considering the density as homogeneous, is justified on the basis of 
the results of Fig. 4, where it is shown that the phenomenon observed in [Q is mainly dependent on the change in 
the density. 

For an ideal gas in a harmonic oscillator potential at T = 0, all the atoms are in the ground state, and a rough 
estimate of the density gives n ~ A^/(47raozaQ^/3). Taking N — 10^ Sodium atoms and Vz — "^Qx 2tt Hz, I'r — 70x 27r 
Hz, rt — 0.567, the ground state dimensions are aoz ~ and apr ~ 2.4/i, and we obtain a density n « 8 x 10^^ 
atoms per cm'^. Thus, according to (55), the group velocity is uj^deal ^ q j^/c^qq^ 



For an interacting gas in the Thomas- Fermi limit, the cloud is an ellipsoid of axes 2RTFr in the radial direction and 
'^Rtfz in the axial direction, where Rtfj is the Thomas-Fermi radius: 



RxFj = \l ;337-2 with j = r,z , (56) 



and is the chemical potential, defined as 



2/5 



A^=^ — ^ , (57) 
2 \ aho 



with as scattering length, Vho = {vr'^zY^^ geometrical average of the oscillator frequencies and auo = y/^hnvho 
corresponding oscillator size. Taking as = 2.75 nm, for the set of parameters of the experiment the Thomas-Fermi 
dimensions of the cloud are Rtfz ~ ^1 Afi and Rtft ~ 13.6/x . Considering the density of atoms as homogeneous, 
we obtain n w 3 x 10^^ atoms per cm^. From Eq. (|5^ ) we find for the group velocity vj"^ ~ 9 m/sec, which is 
comparable with the value measured in |^] for temperatures below the critical temperature. Therefore, for 10^ atoms 
we find a difference of two orders of magnitude in the value of the group velocity between the ideal case and the 
interacting case. Such difference increases or decreases depending on the total number of atoms in the trap. This 
estimate substantiates the inference that interactions are responsible for a lower density, and therefore, for a higher 
average group velocity of the light. 



V. CONCLUSIONS 



We have derived an approximate analytical expression for the group velocity of a pulse propagating through an 
ultracold gas which is confined in a box and by a harmonic potential. We have shown that the results reproduce 
qualitatively the experimental ones presented in 0. From our analysis it emerges that the definite variation of the 
group velocity with the temperature of the gas is an effect related to the variation of the spatial density of the 
gas. We see that the ideal gas model provides a qualitative description of the results for T > Tc- However, the 
behaviour at T < Tc can be described in a satisfactory way only by including the interactions and the fact that the 
cloud is cooled by means of evaporative cooling. The last one has the effect of making the total number of atoms 
temperature-dependent. Such effects will be the subject of future investigations. 
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FIG. 1. Level scheme 



FIG. 2. Plot of the group velocity in m/sec as a function of the effective temperature 6 — T/Tc, for a gas of sodium atoms: 
(a) free bosons (dashed line) and (b) trapped bosons (dotted line). Here Q = O.567, Ag = A,. = 0, V gr = 27rl000 Hz. For the 
free case: Tc = 154 nK, n = 3.8 x lO^^cm"^; For the harmonic potential case: Vr = 1-k y. 70 Hz, Vz — li^ y~ 20 Hz, Tc = 432 nK, 
iV = 8.3 X lO''. 



FIG. 3. Onset: Plot of the group velocity in m/sec in logarithmic scale as a function of the effective temperature Q = T/Tc for 
a gas of sodium atoms as in [7]. The upper curve corresponds to f2 = I.27, whereas the lower curve corresponds to = O.567. 
Here, = 432 nK, iV = 8.3 x 10*^, Pgr = 27r x 1000 Hz, = 27r x 70 Hz, Vz=2-r y. 20 Hz and Ag = A^ = 0. Inset: Plot of 
the low temperature behaviour of the corresponding curves in linear scale. The radius of the section S is R=15 /i. 



FIG. 4. Calculations with two different radii for the section S: R = 15/i (solid lines) and R — ^ KbT /mur (dotted lines). 
The two bottom (top) curves correspond to = O.567 (f2 = I.27). All the other parameters are reported in the caption of 
Fig.3. 
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